Enhancement of bichromatic high-harmonic generation with a high-frequency field 
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Using a high-frequency field superposed to a linearly po- 
larized bichromatic laser field composed by a wave with fre- 
quency u) and a wave with frequency 2u), we show it is possible 
to enhance the intensity of a group of high harmonics in or- 
ders of magnitude. These harmonics have frequencies about 
30% higher than the monochromatic-cutoff frequency, and, 
within the three-step-model framework, correspond to a set of 
electron trajectories for which tunneling ionization is strongly 
suppressed. Particular features in the observed enhancement 
suggest that the high-frequency field provides an additional 
mechanism for the electron to reach the continuum. This in- 
terpretation is supported by a time-frequency analysis of the 
harmonic yield. The additional high frequency field permits 
the control of this group of harmonics leaving all other sets of 
harmonics practically unchanged, which is an advantage over 
schemes involving only bichromatic fields. 
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Within the last few years, the perspective of obtaining 
efficient laser sources in the extreme ultraviolet (XUV) 
regime has led to the proposal of several schemes for con- 
trolling the harmonic spectrum of an atom subject to a 
strong laser field (/ ~ 10 14 W/cm 2 ), using for instance 
additional static fields JjJ , ultrashort pulses |2| , bichro- 
matic driving fields or additional confining poten- 
tials M. These schemes are based on the "three-step 
model" which describes very well the spectral fea- 
tures observed in experiments. These features are the 
"plateau", where harmonics of roughly the same inten- 
sities exist, and the "cutoff", where the harmonic signal 
suddenly decreases in orders of magnitude. According 
to this model, the generation of high-order harmonics in 
strong laser fields corresponds to a dynamical process in 
which an electron leaves an atom at an instant to via tun- 
neling ionization, this electron in the continuum is then 
accelerated by the driving fields, and, if it comes back 
and recombines with the parent ion at a time t\, a high 
energy photon will be generated. The harmonic energy 
is given by £1 = |e | + -Ekin(*o,*i), with £?kin(*o,*i) and 
|eo| being, respectively, the kinetic energy of the electron 
upon return and the ionization potential of its parent 
ion. The cutoff corresponds to the maximum electron 
kinetic energy. By manipulating the different steps of 
the dynamical process, one can in principle control high- 
harmonic generation. 

If the purpose is to increase the cutoff energy, very effi- 
cient schemes exist. They involve either additional static 
fields [Q, or atoms placed in confining potentials ||. In 
both cases one introduces an additional force which mod- 



ifies the propagation of the electron after injection, result- 
ing in an increased cutoff energy. However, such schemes 
may require extremely high static fields p| , or appro- 
priate solid-state materials, whose existence is still under 
investigation |l2] ]. 

Another scheme for manipulating high harmonics in- 
volves bichromatic driving fields By varying the 
frequencies, the relative phase, and the intensities of the 
two driving waves, one can modify the propagation of 
the electron in the continuum, and even the "first step" , 
i.e., the tunneling ionization. In contrast to the previ- 
ously discussed schemes, high-harmonic generation with 
bichromatic driving fields is already experimentally fea- 
sible0§. 

In the bichromatic case, however, no simple expression 
for the cutoff energy exists , so that it is not straight- 
forward to predict whether the plateau can be extended. 
As a general feature, high-harmonic spectra from atoms 
in bichromatic fields exhibit a plateau with a relatively 
complex structure, which may have several cutoffs. These 
cutoffs are given by local maxima of -f?kin(£oj ^1)7 which 
are, however, unequally important for the resulting spec- 
tra. Thus, it may happen that the absolute maximum 
of -Ekinf^O; ti) leads to a much less pronounced decrease 
in the harmonic signal than a local maximum at a lower 
energy position. 

A very illustrative example is a bichromatic field con- 
sisting of a wave with frequency lu and its second har- 
monic [|-§]> The addition of the second driving wave 
causes a splitting in the monochromatic-field cutoff en- 
ergy e max = £o| + 3.17U P , Up being the ponderomo- 
tive energy. As a direct consequence, there is a double- 
plateau structure in the harmonic spectra, with an upper 
and a lower cutoff, whose energies are higher and lower 
than £ max , respectively. The upper cutoff corresponds to 
the absolute maximum for -Ekin^Oj However, it ap- 
pears in the spectrum only as a small shoulder due to the 
relative low harmonic intensity. The lower cutoff, on the 
other hand, is related to a decrease of orders of magnitude 
in the harmonic yield, being therefore the one effectively 
measured in experiments (see, e.g., |5]J^] for a more com- 
plete discussion) . One can explain the intensity difference 
in this double-plateau structure in terms of the width of 
the effective potential barrier through which the electron 
tunnels. This barrier is given by V c s — V(x) — xE(to), 
where E(to) is the field at the electron emission time 
and V(x) the atomic potential. For the upper cutoff, the 
atomic potential is not as much distorted by the field 
as for the lower cutoff. This results in a considerably 
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wider effective potential barrier and, consequently, much 
weaker harmonics. 

In this paper we consider again the case of a linearly 
polarized u — 2uj field, as in However, this time 

our aim is to increase the intensities of the upper-cutoff- 
harmonics close to the intensities of those belonging to 
the lower cutoff, effectively extending the cutoff energy 
beyond e max = ko| + 3.17f p . For this purpose, one must 
provide an additional mechanism for the electronic wave 
packet to reach the continuum, thus compensating the 
weak tunneling ionization. We demonstrate that this can 
be done with an additional driving wave. This third wave 
has a relative high frequency but low intensity compared 
to the bichromatic field. The nature of the third wave 
is such that it alters the electron injection into the con- 
tinuum and, at the same time, it does not appreciably 
modify the ponderomotive energy and the cutoff. This 
field configuration is similar to that used in where 
the influence of the background field on x-ray-atom scat- 
tering processes was investigated. 

We restrict ourselves to a one-dimensional model, 
which still describes high-harmonic generation with lin- 
early polarized fields well in qualitative terms. We solve 
the time-dependent Schrodinger equation 



i Jt \m) = 



L- + V(x)-p-A(t) 



\m) 



(i) 



numerically, for an atom initially in the ground state, 
with binding potential V(x) subject to a laser field 
E(t) = —dA(t)/dt. Atomic units are used throughout. 
The external laser field is taken as 
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Eqx sin(wif) + E i sin(w;t + <j) U ), (2) 



with Eoi, uji and <pu being the field amplitudes, frequen- 
cies and the phases with respect to the first driving wave, 
respectively. In this paper, we choose lo\ = ui, u>2 — 2a>, 
E03 *C Eqx,Eo2 and 3> U). Unless stated otherwise, 
the relative phase 4>is is set to zero. The binding poten- 
tial is taken as 



V G (x) = -aexp{-x 2 /[3 2 ), 



(3) 



which is a widely used expression for modelling short- 
range potentials. The harmonic spectra are calculated 
from the dipole acceleration x — (i/j(t)\ — dV(x)/dx + 
E{t)\ij)(t)) \ 14|. Spectral and time- frequency analysis is 
subsequently done on x (f ) . This latter method has been 
extensively used to extract the main contributions to 
high-harmonic generation within a cycle of the driving 
field |i|Jlq1. These contributions give the electron return 
times corresponding to a particular set of harmonics (see 
P,[l5| for a more complete discussion). 

Time-resolved spectra are computed by performing a 
Fourier transform with a temporally restricted window 
function. For an arbitrary function f(t'), this transform 
is 



■F(f,fi,o-)= / dt'f(t')W(t,t',n,a) 



(4) 



We consider a Gabor transform, which has a gaussian 
window function 

W(t, if, ft, a) = exp[-(f - t'f/a 2 ] exp[ifif'] . (5) 

The usual Fourier transform J-(£l), for which all the tem- 
poral information is lost, is recovered for a — > 00. The 
temporal width a corresponds to a frequency bandwidth 
(To = 2/er. 

The kinetic energy of the electron upon return is taken 

as 



£ kin (f 0) fi) = ^[A(t 1 )-A(t )} 2 . 



The ponderomotive energy is given by 



A 2 {t)dt 
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(6) 
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To first approximation, if the third driving wave is much 
weaker than the others and its frequency is much higher 
than w, the contribution from the additional field to the 
ponderomotive energy ([?]) and to the kinetic energy (j^) 
can be neglected. More specifically, for the parameters 
used in this paper, we observed, after solving the clas- 
sical equations of motion of an electron in a field (j|), 
that for i?03 /E02 < 1 and W3 > 5u the influence of the 
third driving wave on these two latter quantities was not 
significant. 

As a starting point, we shall discuss the existence of 
the enhancement in question. We will restrict ourselves 
to varying the relative phase c/f>i2 between the two low- 
frequency driving waves and the high-frequency field pa- 
rameters. The bichromatic field strengths are similar to 
these in ||, namely Eqi = 0.1 a.u., £^02 = 0.032 a.u, 
which give the intensity ratio Iiujluj = 0.1. For these 
field parameters, the high-harmonic spectrum displays a 
clear double-plateau structure, with a lower and an up- 
per cutoff (c.f. Fig. 1 and Refs. |J|). The lowest fre- 
quency is taken as to = 0.057 a.u., which is typically used 
in experiments. The ground-state energy was taken as 
£o| = 0.57 a.u., which roughly corresponds to the ar- 
gon ionization potential and, unless stated otherwise, we 
took a = 1.15 a.u. and (3—1 a.u. in (||). This gives a 
model-atom with a single bound state. 

The influence of the third driving wave on the har- 
monic spectra is shown in Fig. 1, for relative phases 
012 = and 0i2 = 0.37T and several field strengths i5o3- 
The relative phase ^13 was set to zero in both cases, and 
0J3 was chosen to be ten times u%. Apart from the strong 
enhancement around W3 (near harmonic order 10), the 
main differences between the case with a third wave and 
the bichromatic case (-E03 = 0) are displayed in the so- 
called "upper-cutoff" harmonics, which lie, for the phases 
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in question, roughly at e u = \eq\ + 3.8U P , with the cor- 
responding harmonic orders near N — 62. The inten- 
sity of this group of harmonics can be changed signifi- 
cantly when the field strength of the third wave is in- 
creased. Whereas for <pi2 = these changes are quite 
irregular, for </>i2 = 0.37T they appear as very pronounced 
enhancements, which reach three orders of magnitude. 
They make the intensities of the upper-cutoff and lower- 
cutoff harmonics comparable, effectively extending the 
harmonic production region. This occurs already for a 
third wave whose intensity is only a few percent of that 
of the Dichromatic field. 
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FIG. 1. Harmonic spectra for an atom subject to 
the three-color field (El), with strong- field amplitudes 
Eoi = 0.1 a.u., Eq2 — 0.032 a.u., high-frequency field 
strengths < Eos < 0.02 a.u., frequencies wi = 0.057 a.u., 
u>2 = 0.114 a.u. and 0^3 = 0.57 a.u., and relative phase 
4>i3 — 0, for cj>i2 = (part (a)) and <j>i2 = 0.37T (part (b)). 
The lower cutoff (near harmonic order 50) and the upper cut- 
off (near harmonic order 62) are marked with arrows. The 
lower-cutoff energy depends more sensitively on cpi2 than the 
upper-cutoff energy [5]. 

A particularly interesting feature of the scheme is that 
it affects mainly the group of harmonics in the upper 
cutoff, leaving other groups of harmonics practically un- 
affected. This is in contrast with the purely bichromatic 
case, for which a change in the relative phase (f>\2 leads 
to changes in the whole plateau structure [f|||]. In ana- 
lyzing the enhancement effects, we recall that the pres- 
ence of the high-frequency third wave does not modify 
the ponderomotive energy and the classical trajectories. 



Thus, the three-step model leads us to the conclusion 
that a high-frequency induced process is injecting the 
electron into the continuum. This physical picture also 
explains why the remaining harmonic intensities are not 
influenced by the high-frequency wave. Let us consider 
an electron leaving the atom at the emission time toi, 
so that, at its recombination time tu, the lower-cutoff 
harmonics are generated. At the lower-cutoff emission 
time toi, tunneling is so pronounced that, even in the 
presence of the third driving wave, this mechanism is 
still the dominant path for the electron to reach the con- 
tinuum. Thus, the high-frequency induced process is in 
comparison negligible. A similar argument explains the 
differences between Figs. 1(a) and 1(b). For (f>u = 0, 
there is still enough tunneling to compete with the high- 
frequency induced process, such that the quantum in- 
terference between both processes leads to irregular in- 
tensity variations. For (f>i2 — 0.37T, on the other hand, 
tunneling is strongly suppressed, such that the electron 
reaches the continuum mainly due to the high-frequency 
induced process M. 
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FIG. 2. Time-resolved spectra for the lower (part (a)) 
and upper (part (b))cutoff harmonics, for the same bichro- 
matic-field parameters as in Fig. 1, (f>\2 = 0.3-7T and (j>i3 = 0, 
and several high-frequency field amplitudes Eo3- The thick 
dotted lines correspond to E03 — 0. The temporal width 
was chosen as a — 0.1T, such that the frequency width 
is of roughly three harmonics. The precise center of the win- 
dow function has been chosen as fij = 48w (part (a)) and 
£X U = 58cj (part (b)). The return times tu and ti u are indi- 
cated in the figure. 

This interpretation is supported by the time-frequency 
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analysis. In Fig. 2, we display the time-frequency yield 
for groups of harmonics centered at the frequencies f2; = 
48w (Fig. 2(a)) and tt u = 58w (Fig. 2(b)), for several field 
strengths and a relative phase 4>\2 = 0.3ir. The fre- 
quency Qi corresponds to the lower-cutoff energy and Q u 
was chosen at the most prominent harmonic frequency 
in the upper cutoff. These time-resolved spectra are pe- 
riodic within a cycle T = 2tt/lu of the driving field, ex- 
hibiting peaks at times tu = 1.4T and t\ u = 0.9T. These 
are the electron return times related to the lower- and 
the upper-cutoff harmonics, respectively, and have been 
computed in ||, together with the corresponding emis- 
sion times (c.f. Fig. 1 in there). For a given set of har- 
monics, if the high-frequency field is enhancing the injec- 
tion of the electron wave packet at the emission time to, 
the corresponding peak at t± in the time-resolved spectra 
should get more prominent as Eq 3 is increased. This oc- 
curs only for the upper-cutoff return time t\ u . This can 
be observed in both parts of the figure. 

In the following we understand how the observed en- 
hancements depend on the relative phase 4>\ 3 between the 
lowest- and highest-frequency fields, on the field strength 
E03, on the frequency W3, and whether the atomic poten- 
tial has any influence on it. 

The behavior with respect to 4>i3 gives further infor- 
mation about the nature of the effect due to the third 
wave. In case this wave enhances tunneling ionization 
by changing the effective potential barrier, or increases 
the harmonic intensities by distorting the electron mo- 
tion in the continuum, this influence will be different for 
different phases ^13. We verified that the upper-cutoff 
harmonics remain unchanged when ^13 is varied. This 
result suggests that the third wave is playing a very dif- 
ferent role than the bichromatic field. These results are 
displayed in Fig. 3. 



uration intensity for the upper-cutoff harmonics (in the 
figure, close to E 03 = 0.01 a.u.). Furthermore, these har- 
monics are unequally enhanced, the maximum enhance- 
ment occurring, for the example in question, at the har- 
monic order N — 58. This maximum can be slightly 
displaced, depending on the frequency W3. In Fig. 4, we 
show explicitly the behavior of this maximally enhanced 
harmonic with respect to the field strength -E03, for sev- 
eral frequencies W3. This figure shows that the enhance- 
ment effect increases quickly in the weak field region but 
it eventually saturates when reaching the strong field re- 
gion, which means that upper-cutoff harmonics can not 
be indefinitely enhanced by increasing the strength of the 
high frequency wave. 



E 









co 3 =0.285,N=61 




G> 3 =0.456,N=60 




a 0.570.N 58 




68-1. \ 59 



0.000 0.005 0.010 0.015 0.020 0.025 0.030 
E 03< a - U -) 

FIG. 4. Harmonic yields as functions of the high-frequency 
field strength E03, for the same field as in the previous figures, 
with 0i2 = 0.3vr and 0.285 a.u. < co 3 < 0.684 a.u.. The 
harmonic orders displayed in the figure correspond to the most 
efficient enhancements obtained. 
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FIG. 3. Harmonic spectra for the same bichromatic-field 
parameters as in the previous figures, cj>i2 = O.Stt, field 
strength E03 = 0.015 a.u., and several relative phases (j>i3- 

Another very important issue is the dependence of the 
enhancements on the field strength Eq3- A closer inspec- 
tion of Figs. 1(b) and 2(b) suggests that there is a sat- 
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FIG. 5. Harmonic yields of neighboring upper-cutoff har- 
monics as functions of the frequency ratio U13 /u>i , for a driving 
field as in the previous figures and E03 = 0.015 a.u.. The ion- 
ization threshold is at L03 = lOi^i. 

In Fig. 5 we show the dependence of the enhanced har- 
monics with respect to the frequency W3, for 5wi < LU3 < 
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12u>\. The figure displays major peaks at integer fre- 
quency ratios UI3/U1 and sub-peaks between the integers. 
The enhancements are specially strong at the integers, 
where the frequencies of the three driving waves are com- 
mensurate. The strong enhancements at these particular 
frequencies can be understood as a consequence of the 
high- harmonic generation process: since the production 
of harmonics takes place within many cycles of the driv- 
ing field, one expects it to be most efficient when this 
field is periodic [Q. Fig. 5 also demonstrates that the 
observed enhancements are not related to threshold ef- 
fects, since they occur when the frequency of the third 
wave is either below or above the threshold. 

The investigations presented above are from our nu- 
merical calculations for the short-range potential (^), 
with a and (3 chosen in such a way that it has a single 
bound state. To rule out the possibility that the effect is 
due to an artifact introduced by this model potential, we 
performed similar calculations for several potentials with 
different values for depth a and width j3. We also stud- 

1 /2 

ied the soft-core potential V c = -a [x 2 /(3 2 + l] 1 . We 
observed that the enhancements are present in all situa- 
tions, being however stronger for short-range potentials. 
This is related to the fact that tunneling ionization is 
more efficiently suppressed in the short-range case, such 
that it does not compete with the high-frequency induced 
process. For the gaussian potential used, we obtain very 
pronounced enhancements for (3 < 2 a.u., which is well 
within the experimental range. 

In conclusion, we have shown that, with an additional 
high-frequency field, we are able to enhance a group of 
very high harmonics of a bichromatic driving field con- 
sisting of a wave of frequency u) and its second harmonic. 
These harmonics have energies up to |eo| + 4C/ p , which 
is about 30% higher than the monochromatic cutoff fre- 
quency £ max = ko| + 3.17Z7 P . They correspond to a set 
of electron trajectories for which tunneling is not very 
pronounced in the bichromatic case. The high-frequency 
field provides an additional mechanism for the electron 
to reach the continuum, resulting in the enhancement of 
this group of harmonics. 

This physical interpretation is supported by several 
characteristics of the enhancements, such as, for instance, 
the fact that the remaining groups of harmonics are not 
altered by the high-frequency field. For an electron leav- 
ing at a time to such that the lower-cutoff and plateau 
harmonics are generated, tunneling ionization is the dom- 
inant process. This is true even in the presence of the 
third wave, so that the high-frequency induced process 
plays practically no role. For the upper-cutoff harmonics, 
on the other hand, tunneling is relatively weak, so that 
the high-frequency induced process may compete with it 
or even be far more prominent. This latter effect is ob- 
served in the behavior of the enhancements with respect 
to the relative phase <pi2 between the two low-frequency 
waves. The time-resolved spectra of the lower -and up- 
per cutoff harmonics confirm this physical picture. As 



the high-frequency field is increased, the time-frequency 
yield at the return times t\i corresponding to the lower 
cutoff remains pratically unaltered, whereas the yield at 
the upper-cutoff return time t\ u is enhanced in orders of 
magnitude. Thus, the third field is increasing the electron 
injection in the continuum at the corresponding electron 
emission times to u . Finally, since the enhancements are 
independent of the phase 4>i3 between the highest- and 
lowest frequency fields, they can not be atributted to 
mechanisms such as tunneling ionization or a distortion 
in the electron propagation. 

The enhancements are particularly strong when the 
relative phase between the two strong driving waves is 
chosen such that tunneling ionization is strongly sup- 
pressed. Such a case is provided by taking $12 = 0.3n 
Q , for which an appreciable enhancement is already ob- 
tained with high-frequency fields whose intensities are 
only a few percent of that of the bichromatic fields. This 
phase control is already experimentally feasible An- 
other prerequisite for pronounced enhancement effects 
concerns the frequency of the third wave, which must 
be an integer multiple of the fundamental frequency of 
the bichromatic field. Furthermore, our theoretical stud- 
ies show that the enhancements are always present for 
very different potentials, but they also suggest that a 
short-range potential gives stronger enhancements and 
therefore is a more appropriate choice for a possible ex- 
perimental consideration. 

It is also worth mentioning that, in principle, a simi- 
lar enhancement scheme is applicable to any sets of har- 
monics corresponding to a set of electron trajectories for 
which tunneling ionization is strongly suppressed. For 
the specific uj — 2uj configuration considered in this paper, 
this condition is satisfied by the upper-cutoff harmonics. 
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